Do Binary Hard Disks Exhibit an Ideal Glass Transition? 
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We demonstrate that there is no ideal glass transition in a binary hard-disk mixture by explicitly 
constructing an exponential number of jammed packings with densities spanning the spectrum from 
the accepted "amorphous" glassy state to the phase-separated crystal. Thus the configurational 
entropy cannot be zero for an ideal amorphous glass, presumed distinct from the crystal in numerous 
theoretical and numerical estimates in the literature. This objection parallels our previous critique 
of the idea that there is a most-dense random (close) packing for hard spheres [Torquato et al, Phys. 
Rev. Lett, 84, 2064 (2000)]. 



Understanding the glass transition in dense or super- 
cooled liquids remains one of the challenges of condensed 
matter physics. In particular, considerable effort has 
been directed at identifying the cause of the dramatic 
slowdown of the dynamics in the vicinity of the kinetic 
glass transition, as evidenced in a decrease of the diffu- 
sion coefficient and an increase in relaxation times. One 
possibihty is that a thermodynamic transition different 
from the usual liquid-sohd transition underhes the ki- 
netic one. One scenario originally suggested by Adam 
and Gibbs 01 relates the slow diffusion to a vanishing of 
the number of alternative configurations available to the 
liquid, leading to an ideal thermodynamic glass transition 
when the liquid has no choice but to remain trapped in 
one of few glassy configurations. An important basic as- 
sumption in these considerations is that crystalline con- 
figurations, which are thermodynamically favored, are ki- 
netically inacessible and therefore the liquid is restricted 
to exploring "amorphous" configurations. In particular, 
the term amorphous has become implicitly attached to 
the term glass, and crystalline configurations have been 
assumed to be qualitatively different from glassy ones. 
In this Letter, we study a specific model glass former, 
namely, a binary hard disk mixture, and show that, for 
this model, the presumed "ideal glass" is in fact a phase- 
separated crystal, and that that there is no special amor- 
phous (random) state, but rather a continium of states 
from the most disordered one to the most ordered one 



An inherent-structure formalism was proposed by Still- 
inger and Weber and has since been used extensively 
in the analysis of the thermodynamics of supercooled 
liquids 0- The "inherent-structures" of hard-particle 
systems are in fact (collectively) jammed packings Q, 
which are mechanically stable packings where the parti- 
cles are trapped in a static configuration despite ther- 
mal or external agitation. For soft-particle systems, 
an essential quantity in this thermodynamic analysis is 
the number of distinct energy minima (basins) with a 
given energy per particle. For hard-particle systems 
this becomes the number of distinct jammed packings 
Ngi't'j) = exp [iVsc(</'j)] with jamming packing fraction 
(density) 4>j, where Sc{4>j) is the configurational entropy 
(degeneracy) per particle. It is assumed that the liquid 



remains in the vicinity of these basins for long periods of 
time, jumping from one basin to another as it explores 
the available configuration space. Denser packings are 
favored in terms of their free- volume; the most favored 
one being the crystal of density 0max- However, it is rea- 
sonable to assume that the degeneracy Sc{4>j) decreases 
with increasing (pj. The liquid achieves minimum free en- 
ergy by trading off degeneracy for free volume, so that at 
a given density (j) it predominantly samples glasses with 
jamming density 4>j{(j)). The conjectured ideal glass state 
corresponds to the point where the number of available 
basins becomes subexponential, that is, Sc(0j^) ~ 0. At 
densities above an ideal glass transition density (j)iG, de- 
fined via 0j(0ig) = the liquid becomes permanently 
trapped in the ideal glass state. A crucial unquestioned 



assumption has been that 



< 



i.e., that there is 



a gap in the density of jammed states between the amor- 
phous and crystal ones. We will explicitly show that this 
assumption is flawed for the binary hard-disk mixture we 
study, and suggest that this is the case in other similar 
models, contrary to numerous estimates for cfij^ in the 
literature [IEIjI- 

Previous simulations have cast doubt on the existence 
of ideal glass transitions in hard-particle systems ■ 
In particular, it has already been suggested that the slope 
of Sc{(t>j) at dramatically affects the location of the 
presumed transition; an infinite slope shifts the transi- 
tion to zero temperature 0. Questions have also been 
raised about the validity of extrapolations into tempera- 
ture/density regions that are inacessible to accurate com- 
puter simulations , as well as the impact of finite-size 
effects III this Letter, we present clear evidence that 
the concept of an ideal glass transition is flawed for dis- 
tinctly different reasons. Specifically, for our model, we 
explicitly construct an exponential number of jammed 
packings with jamming densities that span from the 
accepted "amorphous" state with (j)j sa 0.84 to that of 
the crystal, with (j^max = n/^/T2 « 0.91, thus clearly 
showing that the configurational entropy cannot be zero 
for the hypothetical most-dense amorphous (ideal) glass 
distinct from the crystal. This objection is in the same 
spirit as the critique of the concept of random close pack- 
ing (RCP) raised by one of us in Ref. Q, namely, that 
there is a continuous tradeoff between disorder (closely 



2 



linked to degeneracy) and density, making the definition 
of a most-dense random packing ill-defined. Instead, Ref. 
01 replaces RCP with the maximally random jammed 
(MRJ) state, i.e., the most disordered of all jammed 
states. 

Here we study a binary mixture of disks with a third 
(composition xb — 1/3) of the disks having a diame- 
ter /c = 1.4 times larger than the remaining two thirds 
{xA = 2/3). Bidisperse disk packings with this aspect 
ratio and xa — xb = 1/2 have been studied (21| as 
model glass formers For this k, it is believed that 

the high-density phase is a phase-separated crystal. Our 
free-energy calculations predict that at the freezing den- 
sity 4>F ~ 0.775, a crystal of predominantly large parti- 
cles should start precipitating from the liquid mixture, 
i.e., systems in true thermodynamic equilibrium should 
crystallize just like the equivalent monodisperse system. 
Nucleation is however kinetically strongly suppressed due 
to the need for large-scale diffusion of large disks toward 
the nucleus 0|, and in fact, we have not observed crys- 
tallization even in simulations lasting tens of millions of 
collisions per particle well above the estimated freezing 
density. 

The quantity Ng{(j)j) has recently been estimated via 
explicit enumeration for binary mixtures of relatively 
small numbers of hard disks (15j . These studies have ob- 
served an approximately Gaussian Ng{(pj) that is peaked 
at a density 0mrj ~ 0.842, interpreted to correspond to 
the MRJ state for this system. Such a Gaussian Ng{ipj) 
corresponds to an inverted parabola for Sc{4>j), as has 
been assumed in previous studies of the thermodynam- 
ics of binary disk glasses For large systems, such 
enumeration is not yet possible and thermodynamics has 
been used to obtain estimates of Sc{4>j), namely, it is es- 
timated as the difference between the entropy (per par- 
ticle) of the liquid SL{(f>) and the entropy of the "glass" 



SLi't') ~ Sg{(t))- Here sl is obtained 
via thermodynamic integration of the equilibrium liquid 
equation of state (EOS) from the ideal gas limit, while 
Sg is defined as the entropy of the system constrained to 
vibrate around a single basin with jamming density 
without the possibility of particle rearrangements. There 
is significant ambiguity in dcfinining these constraints; 
however, at least in the truly glassy region, the system is 
typically spontaneously constrained (jammed) by virtue 
of a very slow rearrangement dynamics, so that Sg can 
be defined reasonably precisely. 

We have calculated Sg{(f>) using a collision-driven 
molecular dynamics (MD) algorithm that will be de- 
scribed in detail in a future publication For hard 
disks, it is very similar to the tether method of Speedy 
13 . The particles are restricted to remain inside hard- 
wall cells in the vicinity of their initial configuration. The 
cells are initially about twice as large as a particle (in 
linear dimension), which allows particle vibrations but 
restricts particle rearrangements, and has been assumed 
to define Sg. The cells then shrink in size slowly during 
the course of the MD run, and at the end they become 



disjoint, leading to a non-interacting system of particles 
with analytically-known free energy. The work done to 
shrink the cells then gives the initial entropy of the glass. 
Closely related methods have previously been used to cal- 
culate configurational entropy |MI3> with similar, though 
less accurate results. For soft-particle glasses, an alter- 
native method is to use the harmonic approximation to 
the vibrational entropy at an energy minimum as an es- 



timate of Sg 1^ |y| . Both methods are rigorous only in 
the high-density (jamming) or T ^ limit, and are ap- 
proximate for liquids, so the quantitative results at low 
(f> should be interpreted with caution. 

The calculation of the true equilibrium liquid equation 
of state (EOS) is not possible inside the glassy region with 
conventional simulation methods, especially for large sys- 
tem sizes 0, 0, 0| . In our MD algorithm, we produce 
glasses by starting with a low-density liquid and growing 
the particle diameters at a growth rate 7 = dD/dt ^ 1 
01, for a very wide range of compression rates 7. Instead 
of looking directly at the reduced pressure p = PV/NkT, 
we assume that a free-volume EOS holds approxi- 
mately and estimate the jamming density 4ij{4i) from the 
instantaneous pressure 



V = 



l-2/p(0)- 



In the jamming limit {p —* 00), we get that cjjj ^ cf) = (f>j, 
and close to jamming « 0j, and thus it is much 
more useful to plot 4>j{4>) instead of p(^). This is evi- 
dent in Fig. ^ where we show (l>j{4>) for a range of 7's. 
We see that fast compressions fall out of equilibrium at 
lower kinetic glass-transition densities 4>g, and that the 
nonequilibrium glassy EOS is very well described by an 
empirical <pj — {1 + a)4>j — acj), where a ~ 0.133, over a 
wide range oi (p > <j)g. It is also clear that even for the 
slowest compressions « 0.8, so that equilibrating the 
liquid in reasonable time is not possible beyond this "ki- 
netic glass-transition" density. Very long MD runs, with 
as many as 50 million collisions per particle, have failed 
to equilibriate our samples at a fixed (j) = 0.8, and in fact 
very different microstructures all remained stable for very 
long periods of time. 

The final jamming densities of the glasses compressed 
at different rates are shown in Fig. |21 Note that slower 
compressions consistently yield denser packings with no 
hints of the existence of a densest glass. Fast compres- 
sions produce packings that are not truly jammed |j| and 
subsequent relaxation of these systems increases the den- 
sity to around (pj « 0.847. This behavior of our hard-disk 
systems is closely related to the observation that super- 
cooled liquids sample saddle points with the saddle index 
diminishing only below the temperature where even the 
slowest cooling schedules fall out of equilibrium Hi, i.e., 
the kinetic glass transition temperature. Observations 
similar to those in in Fig. [21 have already been made for 
systems of soft particles, e.g., the energy of the lowest 
inherent-structure sampled has been shown to continu- 
ously decrease for slower cooling (3| . 
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Figure 1: The equation of state 4>j{4') for A'^ = 4096 disks 
as observed by compressiuK a liquid with different expansion 
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Figure 2; Final jamming densityt^j for different compression 
rates 7, with and without additional relaxation to ensure a 
truly jammed packing has been reached. 

The measured Sc((/>) for the different glass compres- 
sions are shown in Fig. |3| For comparison, the results 
for a slow compression of a monodisperse system are 
also shown and the entropy of mixing Smix = xa In xa + 
xb Inxs has been subtracted from Sc- It is seen that for 
the monodisperse case Sc — Smix (smix = in this case) 
becomes very nearly zero after the liquid freezes (around 
4> ~ 0.7), indicating a continuous or a very mildly discon- 
tinuous liquid-solid phase transition. More interesting is 
the fact that Sc — Smix also becomes nearly zero for the bi- 
nary glasses around the kinetic phase transition (around 
(j) K, 0.8). This important observation has not been made 
before. It means that the estimated number of packings 
that the liquid samples near the glass transition is very 
close to Smix, which is also the entropy of the uncorrelated 
ensemble of discrete states in which a fraction xa of the 
particles is chosen to be large and the remaining particles 
are chosen to be small. It is interesting to observe that 
the parabolic fit to Sc{4>j) from the work in Ref. [l5| . 
if constrained to equal the mixing entropy at the maxi- 
mum, passes through zero at « 0.9, much higher than 
the extrapolation in and close to the crystal jamming 
density. We note that all measurements of Sc in the lit- 



Figure 3: Estimated Sc(</>) — Smix for monodisperse and bidis- 
perse systems of A'' = 4096 disks, as obtained from (suffi- 
ciently slow) compressions with a range of 7's. 

erature that we are aware of are above or close to Sc near 
the kinetic glass transition, and all estimates of the zero 
crossing of Sc are based on extrapolations beyond this 
point without numerical support 0, • 

Such extrapolations are flawed and, in fact, an expo- 
nential number of amorphous jammed packings exist over 
the whole density range from that accepted as the MRJ 
density 0mrj ~ 0.84 to that of the phase-separated crys- 
tal (^ Tpax ~ 0.91. Lower-density jammed packings also 
exist 15]; however, they do not have thermodynamic sig- 
nificance and thus our simulations do not generate them. 
In our simulations we observe that higher t/ij implies mi- 
crosegregation in the form of increased clustering of the 
large particles. This has been most vividly demonstrated 
in Ref. This observation suggests that one can 

generate denser packings by artificially encouraging clus- 
tering. 

To achieve clustering, we start from a monodisperse 
(k — 1) triangular crystal at pressure p — 100 in which a 
third of the particles has been selected as being "large". 
The large particles then slowly grow in diameter while 
the system is kept in (quasi)equilibrium at a constant 
(isotropic) pressure using a Parinello- Rahman-like varia- 
tion of the MD algorithm When k = 1.4 we stop 
the process and then slowly compress the system to jam- 
ming. By spatially biasing the initial partitioning into 
large and small disks, we can achieve a desired level of 
clustering and higher jamming densities. For this pur- 
pose we use a level-cut of a Gaussian random field (GRF) 
with suitably chosen parameters for a flexible family of 
pair correlation functions . Figure 0] illustrates two 
different jammed packings, one with an uncorrelated ran- 
dom choice of large disks, and another with correlations 
encouraging microsegregation. 

To determine the configurational entropy (degeneracy) 
for a given choice of the GRF parameters, we use a 
recently-developed algorithm for obtaining numerical ap- 
proximations of the entropy (per site) of lattice systems 
j20l |. The algorithm numerically measures the probabili- 
ties of observing a given configuration, which in our case 
is the partitioning of the triangular lattice into large and 
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very close to the mixing entropy. Furthermore, our free- 
energy calculations predict a thermodynamic crystalliza- 
tion well-below the kinetic glass transition, casting addi- 
tional doubt on the search for a thermodynamic origin 
of the glass transition. Although the present study fo- 
cused on the hard-disk binary mixture, the fundamental 
principles are general enough to be applicable to a host 
of related systems, notably, both mono- and bi-disperse 
with hard-core and soft interactions. 
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Figure 5: The measured degeneracy of packings of A'' = 4096 
disks obtained by using different parameters of a random 
Gaussian field with Matern correlations ITg|l . as a function 
of the jamming density. For comparison, we have shown 
Sc{(j> = 0.825) for the three compressions shown in Fig. |3| 

small disks, for small windows (we have used 4x4 win- 
dows) . This can be done by generating sufficiently many 
GRF realizations and counting the number of times a 
given configuration occurs. A Markov expansion is then 
used to approximate the entropy per site Sc- Figure 
shows our results for Sc versus the jamming density 0j, 
for a wide choice of the GRF parameters. The results 
clearly show that in order to increase (f)j one must sacri- 
fice degeneracy (sc)- The figure also shows the first mea- 
sured, rather than extrapolated, estimate of Sc{4>j)- This 
observed Sc{(j)j) only goes to zero for the phase-separated 
crystal state, rather than the hypothetical amorphous 
ideal glass state postulated by extrapolations. 

Continuing on work in Ref. 0|, we explicitly demon- 
strated that the concept of random close packing as the 
most-dense jammed amorphous packing is flawed. By 
trading off degeneracy for density in a continuous man- 
ner, we constructed an exponential number of amorphous 
jammed packings with densities spanning the range from 
the most disordered to most ordered jammed states. We 
explicitly calculated, as opposed to extrapolated, the de- 
generacy entropy for densities well above that of the pos- 
tulated ideal glass transition, and found that the degen- 
eracy is positive for all "amorphous" states and vanishes 
only for the phase-separated crystal. For the maximally 
random jammed state, the found a degeneracy entropy 
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